Quantum fields on closed time-like curves 



J.L. PienaarQ C.R. MyersQ and T.C. Ralplj*] 

School of Mathematics and Physics, The University of Queensland, Brisbane 4072 QLD Australia 

(Dated: December 1, 2011) 

Recently, there has been much interest in the evolution of quantum particles on closed time-like 
curves (CTCs). However, such models typically assume point-like particles with only two degrees of 
freedom - a very questionable assumption given the relativistic setting of the problem. We show that 
it is possible to generalise the Deutsch model of CTCs to fields using the equivalent circuit formalism. 
We give examples for coherent, squeezed and single-photon states interacting with the CTC via a 
beamsplitter. The model is then generalised further to account for the smooth transition to normal 
quantum mechanics as the CTC becomes much smaller than the size of the modes interacting on it. 
In this limit, we find that the system behaves like a standard quantum mechanical feedback loop. 
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I. INTRODUCTION 

The existence in general relativity of solutions that 
contain closed time-like curves has long posed a prob- 
lem to physicists: can the laws of physics accommodate 
time travel without running into paradoxes [TH3]? The 
classic example is the science-fiction scenario in which a 
time-traveller kills his own grandfather, thus preventing 
his own existence and creating a paradox. The problem 
may be stated as 'what happens when we choose initial 
conditions such that the evolution on a CTC contradicts 
those initial conditions' ? The classic solution is an ad- 
hoc restriction on our freedom to choose initial condi- 
tions, known as the 'Novikov consistency condition' or 
colloquially as the 'banana peel mechanism' whereby we 
only allow initial conditions that contain a 'banana peel' 
as a means of ensuring that any would-be grandfather 
killer will slip up and fail at his task[3]. Novikov's idea 
of placing constraints on the initial conditions was origi- 
nally proposed in the context of classical general relativ- 
ity, but it underlies many of the path integral approaches 
to quantum dynamics on CTCs due to Hartle, Politzer 
and others [H |6]. Although these attempts recognised 
that quantum mechanical effects can be significant and 
must be accounted for in CTC space-times, many of them 
relied upon a re-normalisation of the initial state in such a 
way as to exclude paradoxes - a procedure that is reminis- 
cent of Novikov's proposal. As a result, it was found that 
severe problems remain in such theories, arising from the 
fact that the laws of physics in the past are altered due 
to the existence of a CTC in the future j^. The poten- 
tial for new paradoxes is highlighted by the formulation 
of the traditional path integral approaches in terms of 
post-selection (P-CTCs) [7 , for which it has been argued 
that superluminal signalling outside the CTC epoch is a 
consequence . 
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Deutsch[3] was the first to show that quantum me- 
chanics might play a more fundamental role in resolving 
the paradoxes as demonstrated by his toy model of a 
chronology-respecting (CR) qubit (i.e. a point-like two- 
level system) interacting with another qubit trapped on 
a CTC. Instead of relying on a re-normalization of the 
CR-qubits initial state as in other approaches, Deutsch 
drew upon ideas from quantum information theory and 
proposed a consistency condition based on the density 
matrix of the CTC qubit. The resulting model places 
no constraints on the input state, always leads to a 
self-consistent solution (thereby solving the grandfather 
paradox) and does not affect the laws of physics prior 
to the CTC epoch, thereby avoiding many of the prob- 
lems encountered in other models. The fundamental non- 
linearity of Deutsch's solution leads to increased power 
for certain quantum information tasks [101 lllj . Of par- 
ticular relevance to this paper is the re- formulation of the 
Deutsch model in terms of 'equivalent circuits' jl2j. 

The equivalent circuit formulation of the model demon- 
strates that it is possible to keep track of the hidden 
degrees of freedom in the CTC and their correlations 
with the CR-qubit and thereby maintain coherence. In 
the equivalent circuit, the loss of coherence seen in the 
Deutsch model is interpreted as the tracing out of the in- 
accessible degrees of freedom by a detector in the asymp- 
totic future. With this interpretation, the equivalent cir- 
cuit reproduces the results of the Deutsch model. How- 
ever, it also allows us to go beyond the Deutsch model 
in some important ways: first, it explicitly includes the 
CTC degrees of freedom to provide a coherent unitary 
description, which we will use to replace Deutsch's sim- 
ple qubits with more complex field states; second, it's 
structure suggests a further generalisation that allows us 
to treat fields whose wave-packets in space-time are of 
comparable size, or larger than, the CTC itself. In this 
way, we find a smooth transition between non-standard 
and standard quantum mechanics as a function of the 
temporal dislocation produced by the CTC. 

In section II we will review the Deutsch model and its 
equivalent circuit formulation. Section III is concerned 
with the extension of the model from two-dimensional 
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states to those with an infinite dimensional Hilbert space. 
In section IV we discuss a nonhnear modification to quan- 
tum optics motivated by the equivalent circuit, which we 
use to perform the calculations of the previous section in 
situations where the wave-packets become comparable or 
larger than the CTC. Our results are shown to be consis- 
tent with the Deutsch model in one limit and to recover 
standard quantum mechanics in another limit. 



II. QUANTUM CIRCUIT MODELS OF CLOSED 
TIME-LIKE CURVES 

A. Some technical considerations 

Unfortunately, space-times with CTCs do not admit 
foliation into a family of space-like hyper-surfaces, which 
is required for the notion of time evolution of a quantum 
field to be meaningful. Nevertheless, one might imagine 
that the CTCs are confined to a localised epoch in space 
and time, and hope to define some kind of scattering 
matrix between the asymptotic past and future. This is 
the approach used almost exclusively in the literature on 
the topic, and it will be used here. 

As emphasised by Hawking, all such attempts must 
result in a non-unitary scattering matrix if we insist on 
using quantum field theory in its accepted form 1 13] . This 
raises the problem of how to retain a probability interpre- 
tation, since non-unitary evolution seems to imply non- 
conservation of probability. 

It is possible to evade this problem (as Hawking does) 
by pointing out that the loss of coherence comes from in- 
teracting with a part of the universe that is inaccessible 
insofar as it requires a theory of quantum gravity to de- 
scribe its internal dynamics. Since we do not have such a 
theory, we are forced to trace it out, and we are justified 
in doing so as long as we only consider detectors in the 
far future that also cannot access the CTC. The problem 
of whether or not CTCs can be consistently described by 
physics is therefore postponed until we have a theory of 
quantum gravity. A related issue is Hawking's chronol- 
ogy protection coniecture|14|. which provides some indi- 
cation that time travel will be impossible in a final theory; 
however it remains unconfirmed in the absence of such a 
theory. 

It is nevertheless possible to ask: is there any consis- 
tent way to modify the laws of physics to allow time- 
travel without paradoxes? If it can be shown that con- 
sistency of the theory necessarily implies that there is 
no observable time-travel, Hawking's conjecture would 
be confirmed. Alternately, if a consistent model of time- 
travel can be demonstrated through some modification 
of the laws of physics (whose predictions do not contra- 
dict the results of experiments), then the resulting theory 
may tell us something about quantum gravity. 

In the next section we review an apparently consistent 
toy model of time travel, the Deutsch model, which we 
will subsequently generalise and extend to fields. 



B. The Deutsch model 

A qubit described by some density matrix p — 
(assumed to be pure) interacts via some two-qubit uni- 
tary U with a second qubit in the unknown state pcTC 
which emerges from a CTC, as shown in Fig|l] The first 
qubit then enters the CTC, which sends it back in time, 
so that it undergoes the interaction again now playing 
the role of the second qubit. With this interpretation of 
the circuit, the interaction is seen to happen between the 
younger and older versions of the same qubit. 
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FIG. 1. The Deutsch circuit. An ingoing qubit p interacts 
via a unitary U with it's time-traveling partner pcTC- A 
consistency condition is applied to the reduced state pcTC 
entering the CTC in the future at 'F' and emerging in the 
past at 'P'. After solving this consistency condition, a trace 
is performed (signified by the dotted Hue) giving the reduced 
state p' of the qubit at the detector. 

The state pcTC is constrained by the consistency con- 
dition: 

PCTC^TTp,[Um{'iP\®pCTc)U^] (1) 

which can then be used to compute the evolution of j?/;): 

p' = TrcTc[U{\ij){'4>\®PCTc)U^] (2) 

where the final trace is performed over the CTC subsys- 
tem. The resulting input-output map is non-linear and 
non-unitary in general, but always has at least one solu- 
tion for an arbitrary input and a given U . 

The Deutsch model, through the imposition of the con- 
straint [l] amounts to a non-linear modification to the 
laws of quantum mechanics. The nature of this non- 
linearity is different to that encountered in other models 
because of the nature of the constraint. As emphasised in 
|15j , by requiring matching of the density matrix rather 
than matching of the individual states in its pure-state 
decomposition (or equivalently matching of paths around 
the CTC), Deutsch implicitly treats the density matrix 
as an ontologically 'real' object, not just as something 
representative of an observer's state of knowledge (an 
epistemic state). 



C. The equivalent circuit 

The Deutsch model resolves the grandfather paradox 
without placing any constraints on the input state be- 
cause the equation ([l]) always has at least one fixed point. 
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However, for some specific interactions, there is more 
than one fixed point. To choose between them, Deutsch 
was forced to make an additional postulate that singles 
out the solution with the most entropy. Fortunately, it 
is possible to re-formulate the model such that the maxi- 
mum entropy principle emerges as a natural consequence 
of the dynamics. This re-formulation in terms of the 
'equivalent circuit' can be regarded as an extension of 
the Deutsch model that agrees with all of its predictions, 
but has some other nice properties as well. We review it 
briefly here; details can be found in [T^ . 

To obtain the equivalent circuit from Deutsch's circuit, 
we consider the dynamics from the point of view of the 
time-traveling qubit. After passing through the interac- 
tion and entering the CTC on the top rail of the circuit, 
the qubit finds itself on the bottom rail of a new circuit, 
whose top rail is occupied by an identical copy of the 
qubit in its initial state. After the interaction, this copy 
then enters the CTC and goes on to interact with another 
copy, and so on. Formally, we replace Deutsch's circuit, 
FigjTJ with that of Fig[2] 
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FIG. 2. The equivalent circuit corresponding to Fig[T] The 
consistency condition is now enforced by symmetry: the evo- 
lution of a large array of copies of the input state through a 
series of identical unitaries results in a fixed point for the out- 
put. A final trace is performed to obtain the reduced state 
of the outgoing qubit. The evolution up to the final mea- 
surement may be performed in the Heisenberg picture. While 
the predictions are the same as for the circuit of Fig[T] the 
equivalent circuit has a far richer Hilbert space. 

The minimum size of the circuit is dictated by the num- 
ber of iterations necessary to reach a fixed point. The 
existence of such a fixed point is guaranteed by the form 
of the mapf^. This number will be finite when there is 
just a single fixed point; when there are multiple fixed 
points, the introduction of a small amount of decoher- 
ence provides a means of obtaining a single fixed point 
deterministically after a finite number of iterations. Fur- 
thermore, this solution corresponds to the one selected 
by the maximum entropy postulate in Deutsch's model. 
Once the fixed point is found, a final trace is performed 
over the other outputs. The equivalent circuit therefore 
realises all of the predictions of the Deutsch model within 
the framework of a standard quantum circuit, which is 
advantageous because it allows the full quantum tool-box 
to be applied and provides a clear intuition about the be- 



haviour of the CTC for a given choice of unitary. In par- 
ticular, it is clear that the model is free from pathological 
behaviour such as global superluminal signalling, which 
appears in other models jSj [S] . The evolution is still non- 
linear, because of the presence of multiple copies of the 
initial state (in violation of the no-cloning theorem) , and 
it is still non-unitary in general due to the presence of 
the final trace; therefore these intrinsic properties of the 
CTC are preserved by the framework. 

The equivalent circuit also provides a clear intuition for 
extending the model to fields; we turn to this problem in 
the next section. 



III. FIELD OPERATORS 

To generalise the model, let the rails of the equivalent 
circuit represent the mode operators of a field of massless 
scalar bosons (e.g. the quantised electromagnetic field in 
1-D) instead of point-like qubits. In quantum optics, we 
would associate each rail in this circuit with an operator 
of the form: 



A 



G 



dk G(fc, x)ak 



(3) 



representing the annihilation of a photon in the wave- 
packet G{k,x), whose Heisenberg evolution through the 
circuit can be used to define the dynamics. From here 
on without loss of generality we assume the wave-packet 
is a Gaussian superposition of plane waves in flat space: 
G{k, x) — g(k)e*'^^, where gCk) is normalised and is equal 
to zero in the region k < 0. Quantities of interest are 
then obtained by calculating the expectation values of 
the appropriate functions of this operator and its adjoint. 

As we are now explicitly incorporating the space-time 
co-ordinates into the problem, we should clarify what 
metric is being used for the CTC. For simplicity, follow- 
ing the example of Politzer[5], we consider fields evolving 
in flat space-time, where the CTC is implemented by 
making an identification between two space-like hyper- 
surfaces, one of which is in the causal past of the other. 
This is a highly contrived metric in which 'traversal of the 
CTC is treated as a pure temporal dislocation (formally 
defined in terms of translations of the equivalent circuit) 
with no additional dynamics or structure owing to the 
effects of extreme curvature that would normally be as- 
sociated with a 'realistic' metric. The results of section 
III will hold independently of such considerations given 
our assumptions of a large CTC and spatially localised 
modes, which imply that the space-time is locally flat 
along the particle paths. The validity of this argument 
will be re-examined in section IV, in which we seek to 
weaken the assumption of localised modes. 

Returning to the equivalent circuit, each rail is now 
associated with a copy of the original field mode, differing 
only by some index 'm' that delineates the extra degree 
of freedom of the CTC. Formally, we label the rails by m 
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running from — oo to oo. We write: 



A 



G,m 



dk g(k)e''=-ak 



(4) 



for the mode associated to the rrith rail. The modes sat- 
isfy the commutation relation: 



ak. 



(5(k-k')(5„ 



(5) 



in accordance with our requirement that the different 
rails belong to different Hilbert spaces. This leads to 
the same-time wave-packet commutation relation: 



dk.g(k)5*(k) e*'^(''-^')5, 



(6) 



We will find it useful to consider the commutator between 
modes at the same point in space and time, obtained by 
setting X x' in (|6|. We then obtain the 'same-event' 
commutation relation: 



G.m,^O.n 



(7) 



This relation is the main point of departure from stan- 
dard quantum mechanics, because it allows for the pos- 
sibility of interactions between the rails of the equiva- 
lent circuit. Since these rails carry identical copies of 
the input state, interactions between them introduce the 
potential for nonlinear quantum behaviour characteris- 
tic of a CTC. Conversely, in the absence of interactions 
between the rails, the extra degree of freedom becomes 
degenerate and we obtain the commutator of normal, lin- 
ear quantum mechanics independently on each rail; phys- 
ically this is taken to correspond to decoupling from the 
CTC. Our goal is to give a detailed account of the differ- 
ent ways that this decoupling can happen, with reference 
to the particular example of a beamsplitter interaction. 
For the rest of this section we maintain our assumption 
that the dimensions of the wave-packet G(fc, x) are much 
smaller than the CTC; we can then compute the output 
for any reasonable incident field mode, including those 
that contain superpositions of many particles. All we 
have to do is compute the evolution of the modes at the 
detector through the circuit in Fig[2]for a given U, and 
take the expectation values of the desired moments in the 

CO 

initial state = \ipm)- This state represents an in- 

m 

finite number of copies of the original input state, which 
could be any state, such as a Fock state, a coherent state, 
or a squeezed state. We will give examples of all three in 
the remainder of this section, for the case where [/ is a 
beamsplitter. 



A. The beamsplitter on a CTC 

Consider the scenario in which the unitary in FigJT] 
is a beamsplitter; the equivalent circuit for this case is 



shown in Fig|3] The Heisenberg evolution through this 
array is non-trivial because the circuit is formally infinite. 
Fortunately, we can make use of the graph's symmetry 
to compute the evolution by iteration. 
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FIG. 3. (Colour online) The beamsplitter equivalent circuit. 
In the Heisenberg picture, we start with the mode being de- 
tected and evolve it through the unitary operator represent- 
ing the infinite circuit. This produces an expression that is a 
function of the ingoing modes. 



Let the mode at the detector be i 



We wish to derive 



a general recipe for calculating it's evolution through the 
circuit. Following it back along the rail, we see that it 
satisfies: 

where Um+i is the two-mode unitary that acts on the 
(to + l)th and TOth rails. The mode acTCm is given by: 



dCTCm = UrnamUli 
aCTCm-l) 



(9) 

where the last step simply expresses the fact that out- 
put may be written as some function / of the two input 
modes. Iterating this expression, we find: 

aCTCm — /(flm, aCTCm-l) 
= f{am, f{am-l,0,CTCm-2)) 

= /(am, /(am-i,/(am-2, •■•)))• (10) 

Given a specific unitary, we can determine / and hence 
evaluate this expression for dcTCm- Substituting the re- 
sult into ([8| will then give the output mode in terms 
of the input modes. We choose Um to be a linear beam- 
splitter, for which the Heisenberg evolutions for input 
modes a, b are: 



b^ y/Tjb^ 



(11) 



77 a, 



where 77 is the beamsplitter transmittance and </> is an 
arbitrary phase. Then: 



aCTCm-l) 



and ( 10 ) evaluates to: 

00 

acTCm = VvJ2 e"^(v/r^)"t 



(12) 
(13) 



n=0 
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Substituting this into ([8]), we obtain: 

oo 

(14) 

If we set m = and re-label the rails according to 
(ix a-K, which we can do without loss of generality, 
the expression simplifies to: 

00 



n=0 



(15) 



B. The coherent state 



The Heisenberg evolution for a coherent state of am- 
plitude 'a' is: 



= a + a 



(16) 



where b{a) is the unitary displacement operator. The 
initial state is the vacuum: 



|o) = ...|o)i® 10)2 (»...(» |0)c 



(17) 



To calculate the various moments of Oq at the detector, 
we perform the Heisenberg evolution through the equiv- 
alent circuit given by (15) followed by the preparation 



unitary ( 16 ) and then take expectation values in the vac- 
uum state using the commutation relation Q. Since the 
evolution consists entirely of Gaussian operations and 
inputs, we expect the output to also be Gaussian and 
therefore we can characterise the state by its first- and 
second-order momentsjlB]. For the coherent state, we 
find: 



1 - e}<>'^l - rj 



(P)^( — (a't-a')) = l/^/2(7*-7), 



( — (a't+a')) = l/^/2(7*+7), 



Noting that 7 = 6 



Var{P) = Var{Q) 



(18) 



a, together with the results 
( 18 1, shows the output of the equivalent circuit is again a 
coherent state, with the added phase <i>(r/, 0). This is just 
standard quantum mechanics; to a coherent state, a GTG 
looks no more strange than a simple unitary phase shift. 
To explain the intuition behind this result, we recall that 
the beamsplitter interaction between two coherent states 
belongs to a class of problems in which the unitary does 
not produce entanglement between the outputs. This 
means that the function in (|9| has the form: 

f{a.m,acTCm-l) — acTCm-l (19) 



which implies that the interaction can be written as: 

Um — UmUmSWAP (20) 

where acts only on the rrith rail and UmswAP swaps 
the rrith and (m — l)th rails. Substituting this into ([s]) 
gives: 



(21) 



which clearly represents ordinary quantum mechanics 
and is independent of the GTG. Therefore in all situa- 
tions where the beamsplitter does not generate entangle- 
ment between the input mode and the GTG mode, the 
dynamics is decoupled from the GTG and normal quan- 
tum mechanics is restored. We are immediately led to 
wonder what might occur should the beamsplitter gen- 
erate entanglement. One such example is the case of a 
squeezed vacuum input, to which we now turn. 



C. Squeezed vacuum 



The unitary squeezing operator S'(^) is defined as|17j: 



(22) 

where the polar decomposition ^ = re^*^^ conventionally 
defines the squeezing parameter r > and squeezing an- 
gle 6. The squeezing of a vacuum mode a is given by the 
Heisenberg evolution: 



S{S,)aS\£,)^cosh{r)a^ 



-i2e 



sinh{r)aK (23) 



For a squeezed vacuum input, the initial state is again the 
vacuum state and (22) is the preparation unitary. Pro- 



ceeding as before, we evolve the modes at the detector 
through the circuit back to the vacuum and we obtain 
the following results from the first and second order mo- 
ments: 

(a') = 0, 

(P) = (Q) = 0, 
Var{P) = i {cosh{2r) - L(0, d, 77) sinh{2r)) 

Var{Q) = ^ {cosh{2r) + L{(t), 9, 77) sinh{2r)) , (24) 



where 
and 



(1 - 2r])cos{2e) + {r/- lfM{e, cj) 
2 + (77 - 2)77 + 2(77 - l)cos(2(/)) 



M{e, (/)) = cos{2d + 4(f>) - 2cos{2d + 20). (25) 

The function L{<j),9, rj) takes real values between -1 
and 1. If L{(f),9,ri) = cos{29) we recover the result char- 
acteristic of a squeezed state with squeezing parameter r 



FIG. 4. (Colour online). Plot of L{0,e,7j). We see that 
\L{0,Ti)\max = 1 regardless of 77, so there is no added noise. 
Furthermore there is no rotation of the squeezing angle since 
9 is constant with -q. Intuitively this occurs because (j) = 
means that the squeezed states always combine in phase on 
the beamsplitter, leading to no entanglement generation. 




FIG. 5. (Colour online). Plot of L(J,0,77). Some noise is 
present for < 77 < 1, and 6 rotates as a function of rj up to 
+ f at ?7 = 0. 



and squeezing angle 9, i.e. the input state; such a resuh 
would indicate that the state is completely unaffected. 
Unsurprisingly, we find that this occurs trivially when- 
ever 7] — 1, signifying total transmission. In the limit of 
total reflection, ry = 0, we find that L = cos{26 + 20), 
signifying a rotation of the squeezing angle by 0, so again 
the CTC has no effect up to a phase. However, in be- 
tween these two limits, the CTC will not just rotate the 
state but will introduce noise as a function of (j) and 77. 
To characterize the noise, we consider the maximum of 
the absolute value of L, where the maximisation is taken 
over 9 with (j) and rj remaining free parameters; we will 
denote this quantity \L{(f>,9max,r])\ = \Li4',v)\max where 
9max{4'T'n) is the (not necessarily unique) value of 9 that 



FIG. 6. (Colour online). Plot of L(f , 6^, 77). When = | we 
obtain the maximum noise: |I/(f , |)|maa: = 0. This is again 
accompanied by a rotation of ^ by up to = ^. 



maximises |i(0, 0, ry)! for a given and rj. It has the 
property that whenever there is no added noise, there 
will always be some 9max such that 1^(0,77)1^02; = 1- 
In the graphs Figs 4]|6 (colour online) , the presence of 
noise can be determined by noting whether the oscilla- 
tions along the 9 axis are subject to damping; if they are 
not, then there is no noise. 

In the first figure, = 0, and there is no rotation 
of the squeezing angle 0, nor is there any added noise. 
In the last figure, the most extreme case, = 7r/2 and 
there is a rotation of the squeezing angle by an amount 
that ranges from to 7r/2 as 77 decreases from 1 to 0. 
Noise is present for < 77 < 1, with maximum noise 
occurring at 77 = |, in which case the state is symmetric 
and completely thermal (as indicated in Figj6]by the fact 
that |L(0, 7?)|ma2; = there). For intermediate values of 
0, the centre figure shows that the rotation of 9 ranges 
from to 0, with the amount of rotation increasing as 77 
decreases, and the overall noise decreases as decreases, 
becoming zero when also vanishes. 

In summary, the effect of the CTC is to add noise to the 
squeezed vacuum and to change the angle of squeezing. 
Both effects occur whenever ^ n27r,n = 0,1,2... and 
< 7; < 1, otherwise neither effect is seen. 



D. Single photon state 

When the input is a single photon, the usual approach 
would be to place a one-photon Fock state |1) on each rail 

of the equivalent circuit and take expectation values in 

00 

the initial state \ip) = |l)m- This is possible because 

m 

all the results presented in this section can in principle be 
derived in the Schrodinger picture. However, we will find 
it useful for our considerations in section IV to take ex- 



pectation values in the vacuum state ( 17 1 as we have done 
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for the coherent and squeezed states; for this reason we 
adhere to treating the single photon state in the Heisen- 
berg picture right down to the vacuum. The Heisenberg 
evolution of a vacuum mode a to a mode containing ex- 
actly one photon is described by a unitary single-photon 
source, the details of which may be found in [18]. We 
need not reproduce the full expression here, but only take 
note of the following useful properties derived from it for 
a one-photon mode in the vacuum: 



(0|a|0) =0, 
(0|a^a|0) = 1, 
(0|a^a^aa|0) = 0, 
\a, a^l = 1 . 



(26) 



Proceeding as before, we use (15 1 to write the output 



in terms of the input modes and then ( 26 1 to take the 



expectation values in the vacuum state. We then find 
that the output of the CTC has the following moments: 



(a'^a') = 1, 



(27) 



and 

= (a'tfl'ta'a') 
= 4?7(1 — 77) 



(28) 



n,m,p,q—0 



To evaluate this quantity, we note that ( 26 1 implies 



{1 if (n = (?) 7^ (m = p), 
or {n=p) ^{m = q), (29) 
otherwise 

leading to the equivalent expression: 



E (yr^)"+'"+^+^(4aLa,a,) = 

n,m,p,q— 
00 

E (v/r^)"+'"+^+'^ 

n,ra,p,q—Q 

^ {^nq^mp "t" ^np^mq 2 Snm^np^nq} 



This gives us the result: 



J2) 



877(77-1) 



(30) 



(31) 



The result ( 27 1 indicates that the average photon num- 



ber is conserved as expected. For a single-photon state 
however, we would also expect t/^^^ to vanish; while this 
occurs for perfect reflection or transmission, rj = 0, 1, it 
is not the case for val ues of 77 in between these limits. 
In fact, we find from (31) that g'^^ has a maximum of 

24 - 16\/2 sa 1.37... at 77 = « 0.55. This represents 

the value of the reflectivity for which the added noise is 



a maximum. We note that g^^^ has sub-Gaussian statis- 
tics, since the kurtosis of the distribution g^"^^ (77) is found 



to be 



|, which is less than the Gaussian result 



of 3, for all 77. These observations indicate that the sin- 
gle photon state becomes mixed by the CTC, but never 
completely thermalised, unlike the squeezed state. 



IV. GENERALISATION OF THE EQUIVALENT 
CIRCUIT 



In the preceding calculations the particular shape of 
the wave-packet did not play any role, as it was assumed 
to be localised in space to a region much smaller than 
the scale of the CTC. We have seen that in spite of this 
restriction, different choices of input states led to varying 
amounts of coupling with the CTC, and some, like the 
coherent state, did not couple to the CTC at all, collect- 
ing only a phase shift. We would now like to relax our 
initial assumption and ask whether our model can be ex- 
tended to situations in which the modes are longer than 
the CTC, such that the nose of the wave-packet could be 
sent back in time while the tail was still far away. We 
would expect that very long modes might decouple from 
the CTC due to an effective limitation on the interaction 
allowed between the mode and its time-travelling parts. 
In particular, for long modes we expect the noise ob- 
served in section III, in the cases where the unitary was 
entangling, to disappear. In this section, we will modify 
the model of the previous section to take the effects of 
extended wave-packets into account. 

Before we continue, we should address the issue raised 
in section III regarding the role of curvature in our model. 
Thus far we have justified ignoring the effects of curvature 
due to the CTC in two ways: first, we avoid the problem 
of defining modes in a non-globally hyperbolic space by 
treating the formalism strictly as an input-output map 
between asymptotically flat space-times; second, we dis- 
regard the interplay between the spatial properties of the 
modes and the CTC itself by assuming that the modes 
are spatially localised wave-packets much smaller than 
the CTC. It is the latter assumption that we now wish to 
relax, and it could be argued that this calls into question 
the validity of our model. For this reason, the consid- 
erations of this section should be treated as a tentative 
starting point for the construction of a more complete 
model that would include curvature, and the results de- 
rived here as indications of the qualitative behaviour we 
might expect from such a model. In generalising the 
equivalent circuit in this section, we will find it necessary 
to postulate a connection between the spatial and CTC 
degrees of freedom that could point the way for future 
research into the role of curvature in this model. 

Returning to the equivalent circuit, we have not yet 
suggested any physical interpretation for the extra degree 
of freedom. To extend the model, such an interpretation 
will prove useful; let us then consider Figjljin the special 
case where the interaction is the identity. We now find 
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that the incoming mode travels back in time and then 
escapes to the detector without any interaction with its 
younger self. In the equivalent circuit, this scenario in- 
volves nothing more than the detection of an 'older copy' 
of the original mode instead of the mode itself. Thus, an 
incrementation of the parameter 'm' can be interpreted 
as an 'ageing' of the mode. This ageing ought to be quan- 
tified by an invariant parameter depending only on the 
length of the particle's path through the CTC. We there- 
fore select some afhne parameter r that parameterises the 
particle's world line to keep track of the particle's 'age' 
as we have interpreted it; then the incrementation of m 
by some number n corresponds to a shift in the mode of 
uAt along the world-line. This reasoning suggests we re- 
place our model with a more general model in which the 
discrete index 'to' is replaced by a continuous parameter 
ri, whose Fourier complement is t. Then we are led to 
replace Eq|4] with the new wave-packet operator: 



Ag,j= / dk3(k)e^*^ / dr!J(r!,TXn 



(32) 



We will choose the distribution in 51 to be a Gaussian: 
J(ri,r) = j{il)e'^^'^. This choice will be justified from 
physical arguments when we consider the beamsplitter 
example; for the moment we take it as just a mathemat- 
ical convenience. These generalised modes are subject to 
the commutation relation: 



(33) 



We note that a formally equivalent expression to (32) 



was derived by the authors of |19j through similar con- 
siderations; the relevance of that work is discussed at the 
end of this section. Since traversal of the CTC 'n' times 
(or equivalently, a translation of 'n' rails in the equiv- 
alent circuit) produces the transformation J{Q,t) — >■ 
J(ri, T + nAr) in our model, we write down the resulting 
transformed mode as: 



A. 



(n) 



dk5(k)e 



ikx 



dr!j(^^)e^"("+"^")au^n (34) 



It follows from (33 1 that the 'same-event' commutator 
between a mode that has traversed the CTC n times and 
a mode that has traversed the CTC to times is: 



^(n)'^(m) 



dn \j{n)\'e 

(t^-ti)2a 



2 i£l(m— n)Ar 



(35) 



We see that this decays exponentially as the difference 
n — m increases. The rate of decrease of Cn,m is con- 
trolled by the ratio k = ^ where aj is the variance of 

j(r), the Fourier transform of We see that when the 
shift (to — n)AT is much larger than the variance along 
T, the commutator C„_m will vanish when n ^ m and we 
recover ([t]), giving us the equivalent circuit. If, however, 
the commutator C„_m is nonzero for n ^ to, the inter- 
action will be partly decoupled from the CTC. In that 



case we expect behaviour that asymptotes smoothly be- 
tween the equivalent circuit of section III and standard 
quantum optics for which C„.,„ — )• 1. The new mode 
( 32 1 therefore provides the machinery we need to describe 



what happens when the modes become larger than the 
scale of the CTC. To proceed further, we need to estab- 
lish a connection between the function j(t) that deter- 
mines the coupling to the CTC degree of freedom and the 
wave-packet g(x, t) that defines the spatial properties of 
the mode. The nature of this connection becomes clear 
when we consider the beamsplitter example. 



A. The beamsplitter revisited 



It follows from ([34]) and (15) that the output of the 
CTC is now given by: 



n=0 

(36) 

and we replace the 'sharp' commutator ([7| of the equiv- 
alent circuit with the generalised commutator (35). We 



now examine the dependence of our earlier results on the 
parameter k to see what happens when the modes are 
made longer or shorter (in the r direction) compared to 
the CTC. Calculations for the general case are nontrivial, 
but for the limit of ver y lo ng modes, k — ^ 0, we find that 
A(,„) = = A and ( 36 1 becomes the trivial evolution: 



\ n=0 



(37) 



This corresponds to the limit in which the mode does 
not 'see' the CTC due to the variance of j{T) being very 
large. We observe that this is the same result obtained 
in section III in the cases where the CTC became decou- 
pled. However, the derivation leading to ( [37| suggests 
a physical interpretation for the phase shift observed in 
such cases: since the evolution now appears to involve 
just a single mode cycling through the unitary (instead 
of multiple copies of the mode), we note that (37) has 



the same form as a zero-delay feedback loop, because the 
transformed mode at the second input is defined by the 
same operator as the first input mode |20j. To take this 
reasoning further, note that a large feedback loop reduces 
to this same limit when the cross-section along t of the 
spatial wave-packet g(x, i) is very long compared to the 
size of the delay. Hence the limit of very long j (r) in the 
CTC model coincides with the limit of very long ^(t) on 
an ordinary feedback loop with length parameterised by 
t(x, t). The simplest way to account for this coincidence 
is to identify j{T) with g(r) and therefore j{VL) = g{V,) 
for the Fourier transformed modes. This justifies choos- 
ing J(r2, r) to be Gaussian whenever the spatial modes 
are Gaussian. We take this as a postulate to connect 
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the CTC coupling to the spatial properties of the wave- 
packet, leading us to replace (32 1 with: 



Ag,j= / dk5(k)e^'=-y dl]ff(r!)e^"^ak,a. (38) 

The implications of this postulate will be discussed at 
the end of this section; for the moment we merely use it 
as a tool for fixing in our calculations. 



B. Wave-packet decomposition 

For the general case, we expect the output of the CTC 
interaction to lie somewhere between the feedback-loop 
limit and the equivalent circuit of section III. In order to 
perform calculations in the general case, we need to eval- 
uate quantities such as for which the wave- 
packets might only partially overlap. We use the method 
of Rohde, Maurer and Silberhorn (RMS) |21j for decom- 
posing a general wave-packet into components that are ei- 
ther perfectly matched or completely orthogonal to some 
mode of interest. As an example, let the mode of inter- 
est be According to RMS, we can always define a 

complete orthonormal set of functions {^n"*} such that 
An'^ = j4(„) and all the other modes with i ^ are or- 
thogonal to the selected mode. Then we can decompose 
any other mode, say as: 



(39) 



where the operator An contains the accumulated orthog- 
onal modes, and 



A,;= A 



Note that: 



Ao 



1 ^ti 



^(™)'^(n) 



(40) 



(41) 



which is simply the commutator Cn,m of (351 



C. Energy conservation 

Another general result that will prove useful is the ex- 
pectation value {A'^^A'), for any input state. Using (36) 
and separating into matched and orthogonal parts as out- 



lined above, we find: 



X = 



Y = 



jy2 ^ gi(n-m) 
m,n— 



— {n+m) 
T] C„ 



m,n— 



1 - \Cn, 



n=0 



H.c. 



(42) 



The limits of the summations are not easy to determine 
analytically; however, because they converge exponen- 
tially, we can approximate them to arbitrary accuracy 
by truncating after an appropriate number of terms. Af- 
ter doing this, we find that A « 1 and F « 0, leading to 
the result: 



(i'U') = (iU), 



(43) 



regardless of the choice of input state A and independent 
of the overlap Cn,m- This implies that the average num- 
ber of particles is always conserved by the CTC evolution 



(36). Our model of the beamsplitter interaction there- 



fore satisfies global energy conservation for all parameter 
choices, which is an important check of consistency. 



D. Numerical results for different input states 

For the coherent state, as remarked in section III, there 
is no entanglement produced by the beamsplitter and the 
evolution is described by standard quantum mechanics, 
by the application of a phase shift e**'''''^^. We might 
then expect no changes as we smoothly go to the limit 
of a feedback loop, for which k — >■ 0. Performing the 
calculations for different k using the general evolution 
(36), our expectations are confirmed: we obtain the same 



results as (18), independently of the overlap Cn,-, 



For the squeezed vacuum, only the quadrature vari- 
ances differ from the results in ( p4| . This is expected, 
since it is the variances that exhibit the effects of en- 
tanglement and decoherence due to the CTC. Perform- 
ing the numerical calculations, we find that the function 
L((f),9,r]) now depends also on the parameter k, such 
that \L{(j),r])\max ^ 1 when k — > 0, consistent with the 
feedback-loop limit (37). In this limit, the phase shift 



results in a rotation of the squeezing angle as shown in 
FigjTjas a function of ?/, (fi, but there is no noise as there 
is no coupling to the CTC. Plots oi L at cj) = 7r/2 (the 
value of (f) for which the noise is maximised), for k = 1 
and K = 100 are shown in Figs 8]|9 (colour online). As 
K increases, the noise (as measured by the amount of 
damping along 9) increases until we obtain the equiva- 
lent circuit limit. 

Finally, we turn to the single photon state. It follows 
from the result ( [43| that {A'^A') = 1, i.e. the average 
photon number at the output is the same as for the input. 
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FIG. 7. (Colour online). Plot of L{%,e,r)) for «: = 0. We see 
that |-£/(f , ini)\max = 1 for any value of 77, indicating that there 
is no loss of coherence due to the CTC. The overall effect is 
of a phase shift and a corresponding rotation of the squeezing 
angle. 




FIG. 9. (Colour online). Plot of L(f ,6l,ry) for k = 100. The 
decoherence due to the CTC is at its maximum and we recover 
the equivalent circuit limit of section III (compare this graph 
to Fig.([6|). We find lL(f ,r?)U,, = when = |. 




FIG. 8. (Colour online). Plot of L{^,9,v) for k = 1. Now 
l^(f 1 v)\max < 1 for intermediate values of r], signifying some 
thermalisation due to the CTC interaction. 



in this case '1'. The interesting quantity is the probabil- 
ity of detecting photon numbers greater than 1, charac- 
terised by the second order correlation function g'^^\ To 
perform the calculation of f;*-^-* for the equivalent circuit 
in section III, we made use of the identity (30), which 



was derived using the sharp commutato r (|7[ ). As we are 

now using the generalised commutator p5[ ), we need to 

00 ^ ^ ^ ^ 
derive a new identity for the term ^ {Al^Al^ApAq) 

where A is a single-photon mode described by the statis- 
tics (26). Using the RMS decomposition (39), we find: 



n,m,p,ij— 



(44) 



(compare to (30)). Using this result, and truncating the 



summations at an appropriate cutoff, we obtain graphs 
of f;^^-* as a function of rj for the different values of k, 
shown in FigfTO] 




FIG. 10. g*^^' as a function of 77 for the output state when 
the input is a single photon. The graphs represent decreasing 
values of k (light to dark). For k — >■ 00 we obtain the result 
of the equivalent circuit, but as k — >■ the curve flattens out 
to zero, recovering the noiseless feedback-loop limit. 

The graphs shows that for k > 0, there is some proba- 
bility of detecting photon numbers greater than 1 at the 
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output. The input photon might disappear and noth- 
ing come out, or two or more photons might emerge, al- 
though the average photon count must remain 1 accord- 
ing to ([T]). The noise is nonzero only for intermediate 
values of 77, corresponding to the region in which entan- 
glement is created by the beamsplitter (again, there is no 
entanglement and hence no noise when we have perfect 
transmission or reflection). The shape of the distribution 
displays an asymmetry - this can be accounted for by not- 
ing that the physical circuit is itself asymmetric, because 
it takes two reflections for the photon to escape the CTC, 
but only one transmission. Indeed, if we were to alter 
the beamsplitter convention by swapping the outputs, we 
would obtain a mirror-reflection of t;^^^ around 77 — 0.5. 
As K — 0, we approach the result g^^^ = 0, which can 
also be obtained analytically from ( 37 ) . This limit corre- 



sponds to an effective decoupling from the CTC, so the 
noise vanishes and we obtain exactly one photon out with 
certainty. 



E. The role of curvature 

Earlier, we remarked that our generalised model is for- 
mally identical to that found in [T^. There it was con- 
jectured that such a generalised model would be compat- 
ible with CTC interactions - a claim that is confirmed 
by the analysis in this paper. In particular, in that work 
it was suggested that any space-time curvature, not just 
that due to a CTC, should be described by a model of 
the sort that we have introduced here. This would have 
implications for entangled particles in gravitational set- 
tings that could then be tested experimentally; this is 
made possible by the inherent non-linearity of the the- 
ory. It is tempting to disregard such a theory in favour 
of one which reduces to ordinary quantum mechanics in 
the absence of CTCs; however, the physical arguments we 



have made here seem to oppose that view. In particular, 
the equivalence that we have been led to postulate must 
hold between the spatial properties of the wave-packet 
and its interaction with the CTC implies that introduc- 
ing curvature into our model would lead to it becoming 
inextricably linked with the extra degree of freedom; this 
is an interesting avenue for future research. 



V. CONCLUSIONS 

We have seen that it is possible to define a field the- 
ory that is consistent with Deutsch's model for quantum 
evolution on a CTC by applying field modes to the rails 
of the equivalent circuit; this model introduces an extra 
degree of freedom in order to achieve consistency with 
the Deutsch model. We used this model to perform cal- 
culations of the output for a coherent state, a squeezed 
state and a single-photon state interacting with a CTC 
on a beamsplitter. Based on physical considerations, we 
postulated that the extra degree of freedom was related 
to the elapse of an affine parameter along the world-line 
of the time-traveling particle, which led us to a modified 
field theory capable of describing wave-packets compa- 
rable in size to the CTC itself. Using this generalised 
model, we showed that it is possible to smoothly tune 
out the CTC by making the input modes much longer 
than the CTC. In this limit we found that the circuit 
reduces to a feedback loop with zero delay time, as de- 
scribed by normal quantum optics. 
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